Hofstadter butterfly and integer quantum Hall effect in three dimensions 
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For a three-dimensional lattice in magnetic fields we have shown that the hopping along the third 
direction, which normally tends to smear out the Landau quantization gaps, can rather give rise to 
a fractal energy spectram akin to Hofstadter's butterfly when a criterion, found here by mapping 
the problem to two dimensions, is fulfilled by anisotropic (quasi-one-dimensional) systems. In 3D 
the angle of the magnetic field plays the role of the field intensity in 2D, so that the butterfly can 
occur in much smaller fields. The mapping also enables us to calculate the Hall conductivity, in 
terms of the topological invariant in the Kohmoto-Halperin-Wu's formula, where each of (7 xy ,a zx is 
found to be quantized. 



Among the effects of magnetic fields on electronic 
states, one of the most bizzare is Hofstadter's butterfly. 
Namely, when a two-dimensional (2D) periodic system is 
put into a magnetic field, the gap not only appears be- 
tween the Landau levels, but a series of gaps appear in 
a selfsimilar fashion as shown by HofstadterEl. The but- 
terfly refers to an energy spectrum against the magnetic 
flux, <f>, penetrating a unit cell in units of the flux quan- 
tum (f>o = h/e. Usually the butterfly is conceived to be a 
phenomenon specific to 2D. 

Here we raise a question: can we have something like 
Hofstadter's butterfly in spite of, or even because of, a 
three-dimensionality (3D)? This may at first seem quite 
unlikely, since the usual derivation of the butterfly relies 
on the two-dimensionality of the system, so that a mo- 
tion along the third direction (z) should tend to wash 
out the butterfly gaps as well as Landau level gaps. Sev- 
eral authors have.-axtended Hofstadter's problem to 3D 
in the last decadeErtl, and subbands are indeed shown to 
overlap or touch with each other. Here we systematically 
look for the possibility of butterflies (i.e., recursive and 
fractal gaps) in 3D. 

If we do have a butterfly, then we can proceed to ques- 
tion how the integer quantum Hall effect should look like 
on the butterfly. If one examines a theoretical reasoning 
from which the quantization in the Hall conductivity is 
deduced in the usual quantum Hall system, the essential 
ingredient is the presence of a gap in the energy spectrum. 
This was indicated in a gauge argument by LaughlinB, 
and elaborated by Thouless, Kohmoto, Nightingale, and 
den NjisO when a periodic potential exists. There the 
quantized Hall conductivity when the Fermi energy, Ep, 
lies in a butterfly gap is identified to be a topological 
invariant characterizing the position of gaps. 

For 3D Kohmoto, Halperin, and Wu have shown, fol- 
lowing the line of the 2D work, that if there is an en- 
ergy gap in a 3D system, then an integer quantum Hall 



effect should result as long as Ep lies within a gaj 




Montambaux and Kohmoto have actually calculated the 
Hall conductivity in axase where a third-direction hop- 
ping opens some gapsJj If butterflies are realized in 3D, 
we can move on to the systematics of the quantum Hall 
numbers. 

In the present paper we point out that, first, an ana- 
log of Hofstadter's butterfly does indeed exist specific to 
3D under certain criterion that is fulfilled by anisotropic 
(quasi-lD) tight-binding lattices in 3D. In this case the 
butterfly plot refers to energy versus angle of the mag- 
netic field. This is obtained by mapping the 3D system 
to a 2D system. Remarkablly the mapping dictates that 
the ratio of the magnetic fluxes penetrating two facets of 
the unit cell plays the role of the magnetic flux in 2D, so 
that the field intensity B does not have to be strong to 
realize the butterfly. 

More importantly the mapping enables us to systemat- 
ically calculate the Hall conductivity for the 3D butterfly 
via identifying the topological invariant in the Kohmoto- 
Halperin-Wu's formula. We have found that each of 
o~ X y,&zx is quantized in 3D. 

Our model is non-interacting tight-binding electrons 
in a uniform magnetic field B described by the Hamilto- 
nian, 

W = 2(t tfe * 9 «4c i +h.c.) 1 (1) 

(id) 

in standard notations, where the summation is taken over 
nearest-neighbor sites with fy = t x ,t y ,t z along x, y, z, re- 
spectively, 8ij — (e/h) J. A ■ dl is the Peierls phase. We 
first recapitulate the 2D case, because a key in this work 
is a correspondence between 2D and 3D. In 2D with the 
Landau gauge A = (0,Bx), y is a cyclic coordinate, so 
that the wave function becomes ipi m — e Wy7n Fi, where 
(l,m) labels (x, y) coordinates, and v y is the Bloch wave 
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number along y. The Schrodinger equation then takes a 
form of Harper's equation, 

- t x (Fi-x + Fi+i) - 2t y cos(2tt# + v y )F l = EF U (2) 

where <p = Bab/<f>a is the number of flux quanta penetrat- 
ing a unit cell — a x b. While the energy spectrum be- 
comes a butterfly for the ordinary isotropic case, t x = t y , 
the gaps are rapidly smeared out as the anisotropy is in- 
troduced, t y /t x — > 0, since the potential term (the cosine 
function above) weakens. Since t x and t y appear on an 
equal footing, the condition for the appearance of the 
butterfly in 2D is t x w t y . 
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FIG. 1. The energy spectrum of a 3D system with 
t x '■ t y : t z = 1 : 0.1 : 0.1 (a) or a 2D system with 
t x '■ ty : t x = 1 : 0.1 : (b) plotted against the angle 6 in 
a magnetic field (4> y ,(/>z) = 0.2(sin 9, cos 6). 

Harper's equation in 3D can be derived in a similar 
way. For simplicity we consider a simple cubic lattice 
in a magnetic field B = (0, Bsin9, Bcos9) assumed 
to lie in the yz plane. The vector potential is then 
A = (0, Bxcos9, — Bxs'mO), so that y,z are cyclic and 
the wave function becomes ipimn — e Wym+tv * n Fi, where 
(l,m,n) labels (x,y, z). The Schrodinger equation is 

-t x {Fi- X + Fi+i) - [2t y cos(2ir4>J + v y ) 

+2t z cos(-27ny + v z )]Fi = EF h (3) 

where two periodic potentials are now superposed. Here 
(f> y (4> z ) is the number of flux quanta penetrating the side 
of a unit cell (= a x b x c) normal to y(z) (inset of 
Fig.|(a)). 



Although the spectrum in 3D does not in general have 
many gaps (aside from the trivial Bragg-reflection gaps), 
we first notice that a butterfly-like structure does emerge 
for certain choices of (t x , t y ,t z ), as typically displayed in 
|l](a) for (t Xl t y , t z ) — (1, 0.1, 0.1), a quasi-lD system. The 
spectrum is plotted against the angle 9 of a magnetic field 
(4> y ,4>z) = 0.2(sin#, cosd) with b = c assumed here. A 
structure akin to the 2D butterfly is seen in the bottom 
(or at the top) of the whole spectrum. One might con- 
sider this as a 2D butterfly surviving the third-direction 
hopping, but this is wrong as is evident from Fig. 1(b), 
where we turn off t z to find that the spectrum coalesces 
to a series of broadened Landau levels. So we are talking 
about the butterfly specific to 3D rather than a remnant 
of a 2D butterfly. 

We first explore the mechanism why the butterfly 
appears in 3D. The doubly periodic potential in the 
3D Harper equation comprises V^ x \l) oc t y cos{2-K(f) z l + 
v y ), V^ 2 \l) oc t z cos(—2TT(f) y l + v z ). We assume that their 
periods, 1/4> Z and l/<j> y , are much greater than the lattice 
constant (<fi z ,4> y < 1). We also assume that 

ty<\>z > t z (f> y , (4) 

which amounts to an assumption that the local peaks 
and dips of the total potential + is primarily 
that of yW. 

One can then regard the potential minima of as 
'sites', which we call 'wells' to distinguish from the orig- 
inal sites. The wells are separated by l/(f> z and feel the 
slowly- varying V^. Since each well contains many orig- 
inal sites due to the first assumption, we can talk about 
bound states for the well in the effective-mass approach. 
If wells are deep enough, several bound states appear 
and each state forms a tight-binding band (i.e., Landau 
band), and the equation ([|) reduces to 

-tVl'-i + J i'+i) - 2*» cos[-27r(^,/^)Z' 

+ {<py/(f> z )iy y + v z ]J v = EJ V . (5) 

Here t' is a transfer integral between neighboring wells 
labelled by I', Jy the 'effective-mass' wave function, and 
the cosine term is the value of at each minimum of 
V^ 1 ). The reduced equation has exactly the same form 
as that of the 2D system, eqn. (Q) if we translate 

3D :(t x ,t y , t z ,<t> y ,<i> z ) — > 2D :{t', t z , ^ y /<p z ). (6) 

Since the butterfly is a hallmark of an isotropic 2D case, 
one can predict that the subband in 3D for which t' ~ t z 
should exhibit a butterfly. 

We can estimate t' by applying the effective-mass ap- 
proximation to Harper's equation (||). We first convert 
the equation (when there is alone) to a differential 
equation for a continuous variable I = 2-K<j) z l, which turns 
out to contain a combination t y /4> 2 z only (with t x = 1, a 
unit of energy). Since t' is a matrix element of T^ 1 -* cx t y , 
we have a simple scaling law, 
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Note that the value of t' differs from one tight-binding 
band to another, where middle bands, with weaker bind- 
ing, have larger t' . 

We have then numerically calculated t'(t y ,<j> z ) for the 
lowest band. With the scaling, the independence of t' for 
a particular <j) z provides the whole dependence, shown in 
Fig.^j. If we plug in the condition for the butterfly, t z — t', 
the plot indicates how to adjust (j) z to have a butterfly 
for given (t y , t z ). One can immediately find that the but- 
terfly is restricted to the case with t y ,t z <C 1(= t x ), i.e., 
quasi-lD systems. This is because <p z becomes too large 
to satisfy 4> z <§; 1 in the most region of t y ~ 1 or in the 
region t z 1 (out of the plot). 
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FIG. 2. The effective transfer t'(t y , <j> z ) for the lowest Lan- 
dau band in units where t x — 1. By plugging t' —t z , the plot 
serves to identify appropriate values of (f> z to realize the but- 
terfly for given (t y ,t z ). 

This plot shows that the above example, (t x ,t y ,t z ) = 
(1,0.1,0.1), is indeed a right choice, for which we can 
show that t' = 0.05 « t z (= 0.1) for the lowest subband. 
Incidentally, the butterfly is symmetric in this example, 
which is an accident for t y = t z : in Harper's equation 
T^ 1 ) and V 1 -- 2 ^ exchange roles at 6 = 45° for t y = t z . 
Around 45°, or more generally around t y <p z w t z <p y , the 
above argument breaks down but a clear structure re- 
mains. We can explain this as follows. When t y <f) z w 
t z (j) y , exhibits a beat so that the barrier height 

separating the wells varies from place to place, which im- 
plies that t' varies from place to place. However, a change 
in the wall height (ty) changes t' only slightly, since t' has 
a broad peak in Fig.g. 

Thus in 3D the effective flux in eqn.(^) is a ratio <fi y /(j) z . 
In 2D by contrast, a butterfly requires (f> ~ O(l), i.e., B 
has to be impossibly large (~ 10 5 T for a = 2A). Would 
this render the 3D butterfly experimentally feasible?EI 
In principle, Fig.|^ shows that there exists appropriate 
(t y ,t z ) no matter how small <f> z may be. In practice, 
(t y ,t z ) become smaller as <j) z decreases, and the energy 
scale (width of the Landau band 4(i' + t z ) with t' » t z ) 



shrinks for smaller <f> z , so that it will be hard to resolve 
the butterfly structure. For typical quasi-lD organic con- 
ductors such as (TMTSF) 2 X we have t x : t y : t z = 1 : 
0.1 : 0.01 with a,b,c ~ 10A, and we can estimate the 
required <f) z ~ 0.1, which corresponds to B ~ 400T. It is 
still huge, but much smaller than the value required far 
2D and around the border of experimental feasibility.E3 

While we have discussed appropriate values for given 
(t y ,t z ), are there restrictions on (t y ,t z ) to have butter- 
flies? Binding of a well must be so strong that the trans- 
fer to second neighbors is negligible. If one approximates 
a well in with a harmonic potential, the quantized 
energies are (n + \)huj. Then the n-th state is strongly 
bound to each well when this energy is smaller than 4t y , 
the depth of a well. Hence t y should not be too small 
(y/Ty~ > 4> z ), otherwise we have trivial Bragg-reflection 
gaps only. Also, the residual potential V^ 2 ' whose am- 
plitude is 2t z must be weaker than tuv (i.e., t z < ^ zv /t^) 
so that different Landau bands are not mixed. All the 
conditions (those discussed in this paragraph as well as 
fyzii'y ^ Ijty&z 3* t z (f)y) can be interpreted in the 
semiclassical quantization involving the cross sections of 
equipotential surfaces, but the essential hopping (t 1 ) be- 
tween adjacent cross-sectional orbits is outside the scope 
of the semiclassical picture. 

Now we come to our goal of calculating the Hall con- 
ductivity when the Fermi level lies in each gap of the 
3D butterfly. The mapping does indeed enables us to 
accomplish this through identifying the topological in- 
variants in the general formula for the Hall conductivity 
for 3D Bloch electrons by Kohmoto-Halperin-Wu. In the 
formula the Hall conductivity tensor is expressed as 

e 2 x ^ 

°ij = 22 e ij^Gk (8) 

k 

when Fermi energy is in a gap. Here e^-fe is a unit 
antisymmetric tensor, G = nia* + /12b* + /13c* with 
a*,6*,c* being the primitive reciprocal lattice vectors, 
and /xi,^2,M3 are topological invariants specifying each 
gap (i.e., remain constant when we change the direction 
of B in the present context). For an orthogonal lattice we 

have simply <r yz = -^f, <J ZX = ~^f, cr X y = ~Tf- 

So we have only to determine invariant integers, which 
are subject to a Diophantine equation, 

r , P P P 

where we have assumed a rational magnetic flux, 
((j) x ,(j> y ,(f> z ) = ^(n x ,n y ,n z ) (P,Q: integers, n x etc have 
no common divisors), r the number of occupied bands, 
and A another topological invariant. Although the solu- 
tion of the Diophantine equation is not unique, Thouless 
et alu argued for 2D that there is a restriction on the 
integers that decides the solution uniquely. Kohmoto et 
al.u conjecture the uniqueness of the solution in 3D in 
analogy with the 2D case, where the restriction for 3D 
reads \n\n x + ^2% + M3^| < Q/2. 
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FIG. 3. (a) The conductivity (a xy ,a zx ) = —e 2 /h(n3,H2) 
is plotted on a 3D butterfly, where we display the topological 
invariants (^,3,^2) for each gap. (b) The corresponding plot 
for 02D = —(e 2 /h)t on the 2D butterfly. The area enclosed 
by a dashed line in (a) corresponds to the 2D butterfly. 02 in 
(a) corresponds to t in (b), while fj,3 in (a) to s in Eq.rffOl) 

We can then calculate the Hall conductivity for the 
3D butterfly. We assume (^) x ,<p y ,(f) z ) = P/Q(0,n y ^n z ) 
and t y (j) z S> t z tfi y , for which the effective flux in eq. (J5|) is 
4> = 't'yl'Pz = n y /n z . Hence each Landau band should 
split into n z butterfly subbands. Let us consider the sit- 
uation where Ep lies just above the mth subband in the 
lih Landau band from the bottom, i.e., (ln z + m) sub- 
bands altogether. Each subband is shown to comprise P 
bands, so that the gap has an index r = (ln z + m)P. 



Substituting this in the Diophantine eq.(||), we have 
(ln z + m)P = QX + Pn y H2 + P n z^3- Since P and Q 
have no common divisors, s must be a multiple of P, and 
with the above restriction one has A = 0, and we end 
up with ln z + m — n y [i2 + n z /J-3 for the 3D butterfly in 
the lower half of the entire band. Thus we can deter- 
mine /Lt2,M3 for an arbitrary gap in the 3D butterfly as 
explicitly displayed in FigJ|(a). 

If we compare with a corresponding plot for 2D in Fig- 
ure [|(b) , we recognize a beautiful consequence of the 2D- 
3D mapping established here as a unsuspected one-to-one 
correspondence between the Hall conductivities on 2D 
and 3D butterflies as a whole (i.e., for a set of topological 
invariants attached to the recursive gaps). Namely, the 
Hall conductivity in 2D| is given by <t 2 d = — 77 1, where 
t is an integer in a 2D Diophantine equation r = qs +pt. 
If we compare this with the 3D Diophantine equation, 
m = n z (fj,3 — I) + n y /j,2, the mapping dictates a corre- 
spondence n y <~+ p, n z <-> q, m <-> r, so that the invariant 
integers should translate as 

^ 3 - I < — > s, [i 2 < — > t. (10) 

This implies that o~ zx in 3D plays the role of o~2D- 

Future problems are to extend the effective theory to 
the the arbitrary orientation of B, for which Harper's 
equation in the [arbitrary B has been discussed in an ex- 
isting literatureoa. We wish to thank Mahito Kohmoto 
for discussions. 
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